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Summary 

The presentation of aerial radiation patterns in the form of contour plots is 
described and a brief account is given of the plotting algorithm. 
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1. Introduction 

It has been common practice in the past to represent an 
aerial radiation pattern in terms of separate horizontal and 
vertical patterns, it being tacitly assumed that the field in 
any other direction would be given by the product of the 
two. A common form of presentation is the polar plot on 
which field strength at a constant distance is represented 
linearly on the radius co-ordinate. An example is shown in 
Fig. 1 of a horizontal pattern of this type. The polar/ 
voltage diagram is not very satisfactory when it is desired to 
show the performance of a very directional aerial. In Fig. 1 
it will be seen that the field in the minima cannot be read 
with any great accuracy and to overcome this difficulty the 
polar plot is sometimes shown with a logarithmic radial 
scale. A better arrangement is to use cartesian rather than 
polar co-ordinates for the logarithmic plot. This has the 
advantage that the angular scale can be chosen to suit the 
aerial; in particular an expanded scale can be used for an 
aerial with a very narrow beam. Fig. 2 shows the radiation 
pattern of Fig. 1 in this form. This diagram is used in 
connection with the automatic measurement of aerial 
radiation patterns. 
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Fig. 1 - Polar-voltage diagram (u.h.f. panel aerial) 
(RA-114) 



The use of two sections at right-angles (i.e. the radia- 
tion patterns in the horizontal and vertical planes) to 
represent a three-dimensional radiation pattern is acceptable 
only in special cases. With a complex array having a multi- 
lobed radiation pattern a different presentation must be 
used. A diagram which is sometimes used to illustrate the 
performance of large high-frequency arrays is the power 
distribution diagram (sometimes called an 'onion diagram') 
shown in Fig. 3. The characteristic shape arises because the 
length of the lines of latitude are proportional to the sine of 
the latitude; as a result, areas on this chart are equal to the 
corresponding areas on the surface of the corresponding 
sphere. The radiation pattern of the aerial is shown in the 
form of contours of equal field strength or power flow; if 
the latter are plotted the diagram may be used to deduce 
the gain of the array by graphical integration. 



2. Radiation patterns of UHF aerials 

Certain difficulties arise in the presentation of radiation 
patterns of u.h.f. aerials. First, the aerials are often 
several wavelengths long so that there is considerable detail 
in the patterns. Second, the maximum field in any vertical 
plane is normally arranged to fall at some arbitrary angle 
below the horizontal, known as the beam-tilt angle. Where 
a 'horizontal' radiation pattern is required it is usually 
measured or calculated at the beam-tilt angle. However, at 
some stations the beam tilt angle may vary with azimuth so 
that an adequate representation of the three-dimensional 
radiation pattern would need several 'horizontal' radiation 
patterns corresponding to several angles in the vertical 
plane. Diagrams comprising several such patterns super- 
posed have been produced for particular aerials but they are 
not easy to interpret. It was this consideration together 
with the need to be more specific about the detailed 
structure of radiation patterns which led to the devising of 
a new form of representation. 



3. The polar-contour plot 

A number of basic requirements were laid down for the 
new form of presentation. First, all the relevant infor- 
mation relating to the radiation pattern had to be condensed 
onto one diagram. Second, the diagram needed to have a 
pictorial quality so that the reader could readily visualise 
what was represented. Third, the diagram needed to yield 
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Fig. 2- Cartesian-logarithmic diagram (u.h.f. panel aerial) 



decibel values of field strength readily so that it would be 
easy to use in predicting field strength. Finally, it was 
desired to be able to produce the diagram directly from a 
computer. This latter requirement meant that diagrams 
could be considered which would have been unduly time- 
consuming to produce by hand. 

The form of diagram* meeting these requirements is 
shown in Fig. 4. The polar angle corresponds to the angle 
of azimuth with the cardinal points of the compass marked. 
The radial distance corresponds to the angle below the hori- 
zontal, the outermost circle corresponding to the horizontal. 
The range of vertical angles shown (0° to 20°) is chosen to 
include those angles of greatest interest at the majority of 
stations. 

Fig. 5 shows the representation of a typical aerial. 
The contours are of equal field strength and are shown in 
decibels below the maximum. The whole diagram may be 
considered as representative of the way the aerial distributes 
energy over the service area. 

A more complicated example is given in Fig. 6; the 
diagram shows clearly areas of low field strength which 
might easily have been overlooked with other forms of 
presentation. 



First proposed by R.I. Black. 



4. The plotting algorithm 

It was required to produce a plot of contours of con- 
stant effective radiated power from a matrix containing 
samples of the radiation pattern for the appropriate range 
of azimuthal and zenithal angles. Such a programme could 
be adapted to plot contours of any function of the form 
Z = f(d, 0) where Z is the 'height' of a contour. 

The task of devising a contour-plotting algorithm is 
full of pitfalls for the unwary. The basic problem is to 
ensure that no contours of a given height are overlooked in 
the area under construction, that no contours are drawn 
more than once, and that those conditions are fulfilled 
without consuming unacceptable amounts of computer 
storage and run- time. 

Consider a rectangular area ABCD (Fig. 7) divided up 
into a rectangular grid. It is assumed that the values of Z 
are known at every intersection-point on the grid and that it 
is required to plot all the contours of a 'height' Z 
the area ABCD. The procedure is as follows: 



Zq over 



i) Mark every intersection where Z>Z^ with the label 
'high'. Mark all others 'low'. 

ii) Scan every horizontal side of the elemental rectangles 
(such as ab and cd) except those on the boundaries 
AB and CD. Whenever a 'high' intersection has a 



BBC RD 1973/18 



BBC RD 1973/18 




I 

I 



I 

."J 






.■a 



I 




Contour numbers are dB 
below an e.r.p. of 



BBC RD 1973/18 



Fig. 4- Graticule for polar/contour plot 
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Contour numbers are dB 
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Fig. 5 - Polar/contour plot: 4 tiers eachi of two panels set at 120°, 2-5° beam tilt 
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Contour numbers are dB 
below an e.r.p. of lOOwotts 



103.31 BBC UHF PANEL flERIflL CHA1 



Fig. 6 - Polar/contour plot: 4 tiers each of two panels on and 90 ETN. Panels on carry unit power 
and have beam tilt of 2°. Panels on 90° carry quarter power and have beam tilt of 12° 
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F/g'. 7- Method of drawing contours 

'low' intersection immediately to the left of it, mark 
tlie 'higli' intersection as 'uncrossed'. 

Contours starting and finishing on the boundary 
ABCD are now drawn, as follows. Examine the inter- 
sections of the grid-lines with the boundary ABCD 
until there is found a pair for which, from the view- 
point of an observer outside the boundary, the right- 
hand intersection is marked 'high' and the left-hand 
intersection is marked 'low'. P^ and Pj in Fig. 1 are 
such a pair. The start of the contour is obtained by 
inverse interpolation of Z^ between the values of Z at 
P and P . If P4 is 'low' the contour is assumed to 
exit P^ P3 P4 Pg via the side P3 P^. If not, then P^ Is 
examined to see whether the exit is via P^ P^ or P^ 



On whatever side the exit occurs, a new P^ P3 P,, Pj 
rectangle is set up and the process repeated until the 
boundary ABCD is reached. 

In tracing out a contour, whenever a horizontal edge 
is traversed upwards, then the label 'uncrossed' at the 
right-hand end is changed to 'crossed'. The remainder 
of the boundary ABCD is scanned systematically for 
starts of other contours. 

Having drawn all the contours which terminate on the 
boundary ABCD ('open' contours), it is now necessary 
to find those which do not ('closed' contours). This 
is done by scanning the internal horizontal edges 
(such as ab or cd) until one is found for which the 
right-hand end is marked 'uncrossed'. A test rectangle 
P^PgP^Pg is then erected on top of this edge and the 
contour followed in the same way as for open con- 
tours. However, the test for terminating a closed 
contour is that a horizontal edge is encountered such 
that its right-hand end is marked 'crossed'. The rest 
of the internal horizontal edges are scanned and closed 
contours plotted until there are no 'uncrossed' edges 
left. 

The whole of the above procedure is repeated for 
each different contour height. 



Note that the contours are shown as being made up of 
straight lines joining the intersections of the contours with 
the grid. This has the advantages of simplicity and the 
assurance that contours will never intersect each other 
(although they may touch each other). "^ Judicious choice 
of cell-size in the grid makes the straight line method satis- 
factory for most contours. It is thought that the expense 
of curve-drawing procedures would not be justified. 



5. Calculation of the radiation pattern matrix 

The radiation pattern is evaluated from the formula 

i 

where l^".-! is the field intensity in direction (O' 0j') 

(0.', 0j') are a set of directions in spherical-polar 

co-ordinates 
(Ij, nrij, nj) are the direction cosines of (d^, 0j') 
ylje'^i is the current in the ith source 
i?lje'^ii is the radiated field intensity of source 

no. i in direction (9j , 0j ) 
(X|, jj, Z|) are the cartesian co-ordinates of the posi- 
tion of source no. i. 

It is assumed that for each source the co-ordinate 
(X|, y-y Z|), the current ^|e'°i, the direction of fire (0|, (j>{) 
are given and the radiation pattern /?e'^ is given (relative to 
the 'natural axes' of the source) as a function of (0,0). If 
the source is 'turned' through an angle 0| and 'tijted' 
through an angle d-^ (see Fig. 8), then to determine i?e' in 
direction (6-' 0j') it is necessary to evaluate the co-ordinate 
(0|j 0jj) of the direction (0j', 0j') relative to (01,0,). The 
appropriate formulae were determined using the formulae 
of spherical trigonometry and are given below: 

cos 0|j cos 0|j = sin 9j' sin S, + cos 0j' cos 0-, sin (0j' - 0,) 

cos 0|j sin 0jj = cos B^ sin (0j' - 0,) 



sin e^ 



in Q-' cos 0| - cos 0j' sin 0| cos(0J - 0, 



sin 



For the purposes of this programme the axes and co- 
ordinates are slightly non-standard. OX, OY, OZ are 
mutually orthogonal cartesian axes with OZ pointing verti- 
cally upwards. The spherical polar co-ordinates (6 0) are 
such that is measured from OX towards OY and lies 
between 0° and 360°, and 6 is measured from the horizontal 
towards OZ and lies between -90° and -1-90°; +ve angles 
are above the horizontal and -ve angles below. 

It should be noted that individual sources could have 
a further angular degree of freedom apart from 'turn', 0, 
and 'tilt', G. This might be described as 'twist' about the 
direction of fire. This parameter has been excluded from 
the present programme because it would be necessary to 
evaluate both planes of polarisation. 

The most commonly used source at low power relay 
stations is the printed panel aerial.'' Radiation pattern data 
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Fig. 9 - Conversion of polar to rectangular grid 



X Y Z — cartesian axes 

OX'Y'Z' -"natural axes" of source 

{d,<j>) - direction of fire of source 

OP{8',4>)— direction in wtiich RP is to be calculated 

(9"c^") of OP relative to OX'V'Z' is to be determined 
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Fig. 8- Calcuiation of relative 6, 

for this aerial at a number of frequencies are stored in the 
computer along with the programme. This results in a 
great reduction of the number of data cards to be prepared 
for each calculation on arrays of this aerial. 

6. tmplementation of plotting algorithm 

The algorithm described in Section 4 has been followed 
broadly in the contour-plotting subroutine but there are 
some minor differences. In particular, the area over which 
contours are drawn is circular, not rectangular. For the 
contouring algorithm, the polar pattern is cut along the 
North radius and opened out into a rectangle as shown in 
Fig. 9. 

Another difference is that, when the grid is being 
scanned for contour starts, each part of the grid is examined 
for each contour height in turn before moving to the next 
part of the grid. As a consequence, marks 'high', 'low', 
'crossed' and 'uncrossed' must be stored for each contour 
height. 



Provision is made for plotting two sizes of the dia- 
gram. The smaller is contained within A4 size paper and is 
intended for record purposes. The larger, which is almost 
twice the size, is intended for those cases where detailed 
work has to be done using the diagram. Either diagram, or 
both, may be plotted as required. 

The programme is written in Fortran V.* 

Figs. 5 and 6 show the most commonly used ranges 
of vertical and horizontal angles. Any other vertical range 
extending over not more than 20° may be provided and the 
azimuth angle may be restricted if so desired. The cell 
size of the grid must be chosen sufficiently fine in relation 
to the expected complexity of the radiation pattern to 
avoid overlooking significant detail. Intervals of 1° in the 
zenithal angle and 5° in the azimuthal angle were chosen 
for Figs. 5 and 6. In practice, a plot in which the straight 
segments of the contours are not unduly obtrusive is un- 
likely to omit significant detail. 

The computer programme is arranged to draw the 
graticule in a different colour to the contours. If, however, 
it is desired to produce many copies of a diagram, it is more 
economical to produce a plot of the contours only and then 
to print this onto blank graticules. In order to facilitate 
this operation, the letters labelling the cardinal points of 
the compass are drawn by both the contour- and the axes- 
plotting routines and are then available as registration marks. 
Figs. 5 and 6 have been reproduced in this way. 



7. Conclusions 

A new form of aerial radiation pattern diagram has been 
devised to give more complete and accurate information of 
aerial performance. A computer programme has been 
written in Fortran V which calculates the aerial radiation 
pattern and plots it in the new form. 



*The programme calculating ttie radiation pattern matrix and the 
first version of the plotting programme were written by R.E. 
Davies. 
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